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A quantum mechanical model based on a Green's function approach has been used to calculate the 
transmission probability of electrons traversing a two-dimensional electron gas injected and detected 
via mode-selective quantum point contacts. Two-dimensional scattering potentials, back-scattering, 
and temperature effects were included in order to compare the calculated results with experimentally 
observed interference patterns. The results yield detailed information about the distribution, size, 
and the energetic height of the scattering potentials. 
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I. INTRODUCTION 

The electron mobility in a two-dimensional electron 
gas (2DEG) is limited by electron-impurity scattering 
jl], H . The impurity potential originates mainly from 
ionized donors located in the remote doping layer. How- 
ever, the observed high mobilities can only be explained 
when correlations of different donor states are taken 
into account ^, |). Up to now experiments are rare 
which probe the actual impurity potential microscop- 
ically. Standard transport experiments yield only in- 
formation about the averaged potential fluctuations by 
defining the average electron mean free path. Recently, 
it has been shown that the electrons injected via quantum 
point contacts (QPC) into a 2DEG can be used as a local 
probe for these scattering potentials either through the 
interference patterns observed in the transmission proba- 
bility H or through structures in the position dependent 
electron reflection J7I |) . 

In this paper we use a Green's function approach to cal- 
culate the transmission probability of electrons injected 
via QPC into a 2DEG exposed to a small magnetic field. 
The results are compared with experimentally observed 
transmission probabilities where the electrons were de- 
tected with a second QPC in a distance of 4 /jm opposite 
to the injector QPC. In order to reproduce the experi- 
mental results the shape, size, and height of scattering 
potentials has been modeled and back-scattering as well 
as temperature effects have been included into these cal- 
culations. 

II. MODEL 

Our model is adopted to the experimental situation 
presented in Fig. I a) (c.f. Ref. |6)). This figure shows 
two opposite QPCs in a distance L from each other. In 
the experiment these QPCs are formed electrostatically 
in the 2DEG of GaAs/Al .33Gao.67As-heterostructure by 
externally controlled Schottky gates (gray areas). Due to 
the saddle-point-like electrostatic potential of the QPCs 
in the plane of the 2DEG the electrons are injected into 
the region between the two QPCs in form of a collimated 
beam J9|. A weak magnetic field, applied perpendicular 



to the plane of the 2DEG, deflects the electrons and the 
resulting signal of the electrons which reach the detec- 
tor QPC is proportional to the profile of the propagated 
electron beam. It has been shown in Ref. |g) that the 
observed structure in the electron beam signal is due to 
electron interference effects, which originate from scat- 
tering at potential fluctuations imposed into the 2DEG 
from donors in different charged states. A typical exper- 
imental result is shown in Fig. || (curves E). 

The first attempt to calculate the transmission proba- 
bility of an electron beam injected by a QPC into a 2DEG 
and detected by a second opposite QPC was made by 
M. Saito et al. jIo| using a Green's function approach. 
This approach was previously extended by us in order to 
include interference effects due to electron impurity scat- 
tering The scattering potential has been treated in 
a very simplified way by setting the part of a wave func- 
tion at the impurity position to zero (infinite impurity 
potential), which explained the experimental data qual- 
itatively. Here, we present quantum mechanical calcula- 
tions which include a realistic two-dimensional impurity 
potential model as well as back-scattering and temper- 
ature effects. A fitting to the experimental data yields 
information about the distribution, size and strength of 
scattering potentials in a 2DEG. In contrast to our pre- 
vious work Q , finite potential fluctuations with positive 
and negative heights (counted from the bottom of the 
conduction band) are considered, corresponding to the 
regions of reduced and increased electron density, respec- 
tively. Moreover, the back-scattering effects from the im- 
purity potentials and sample boundaries were taken into 
account. It turns out that the back-scattering affects the 
interference patterns considerably in two cases: first, if 
the impurity is located close to the injector or the detec- 
tor QPC (within the phase coherence length), or second, 
if two impurities are located at a short distance from each 
other. The effect is stronger when the scattering poten- 
tial is significantly higher than the Fermi energy. Finally, 
temperature effects are considered, i.e., a thermal aver- 
aging of the propagated beam is included (llj , as well as 
the contribution from electron-electron scattering [[l2| . 

The following initial parameters are set according to 
the experimental situation for the calculations: the dis- 
tance (L) between the two QPCs is 4 /im, which is smaller 
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FIG. 1: a) Scheme of the sample structure used for the calcu- 
lations (^o - the wave function at the exit of injector quantum 
point contact); b) the shape of the scattering potential. 



than the elastic mean free path, the electron mobility and 
carrier density are /i 100 m 2 /Vs and n s ~ 2x 10 15 m~ 3 , 
respectively, and the width (W) of the QPC's exit is 
100 nm accounting for a single mode electron injection 
in case of an adiabatic expansion §,§. The wave func- 
tion propagation is perturbated by potential fluctuations, 
which are described in the following way: Xi, yi, and Axi, 
Ay.i define the position and extension in the x-y plane 
and V(x,y) the potential height (Fig. [I]). It was found 
that in most cases circular impurities (Axj = Ay.;) are 
sufficient to produce a good agreement between the cal- 
culations and the experiments. Hence, we restrict ourself 
in this paper to the presentation of results considering 
only circular impurity potentials with a radius which is 
given by An = ^Ax 2 /A + Ay 2 / 4. 

Our model numerically solves the time-independent 
Schrodinger equation for a static magnetic field [|l0| . Us- 
ing the Green's function method with Dirichlct's bound- 
ary conditions it is possible to evaluate the wave function 
for a line at any distance x' from the injector, according 
to [HI: 
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Here, is the wave function at the injector QPC exit 
which includes the effect of the magnetic field in the form 
of the phase shift Oo between the center and any point 
on the QPC exit: 



* B (0,y) = tfo(0,i/)exp[ie (i/)], 



(2) 



'Fo is the wave function at the injector QPC exit with- 
out magnetic field and G + is the Green's function calcu- 
lated in a weak-magnetic-field approximation using the 
mirror-image method. For a detailed description of this 
method we refer to Ref. p0[. The unperturbated wave 
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FIG. 2: Measured electron beam profiles (curves E) and cal- 
culated interference patterns (curves S) for two scattering po- 
tentials versus magnetic field (simulation results are presented 
with a shift). Parameters for the scattering potentials are: 
a) Xi = 1.3 /im, m — —0.34 /tm, Ar, = 0.045 /im, Vo w 
33 meV for first impurity and Xi = 1.3 fim, yi = 0.18 ^m, 
An = 0.075 /im, Vo ~ 16 meV for second; b) X, = 2 /im, 
j/i = —0.42 fim, Ari — 0.05 ^m, Vo fs 8 meV for first im- 
purity and Xi = 2 fim, yt — 0.215 fim, Ari = 0.075 fim, 
Vo ~ 10 meV for second. 



function can be calculated for the whole 2DEG area 
between injector and detector. However, in the region 
Xi - Ari < x < Xi + Ar t , y { - Ar, < y < y t + Ar t the 
wave function propagation is perturbated by a scattering 
potential. One suitable approximation for the shape of 
this potential is a hyperbolic function |L3| , which is very 
similar to a Gaussian potential commonly used in liter- 
ature for the description of the effective potential due 
to remote charged donors in the 2DEG region 
hyperbolic potential has the advantage that the solution 
for the transmission and reflection coefficients can be ob- 
tained analytically |l3| . The scattering potential is given 
by [Fig. | b)]: 



V(x, y) — Vq cosh 



Viv - Vi) 2 + i x - Xi) 7 
aAri 



(3) 



Vo defines the height of the potential which can be pos- 
itive and negative depending on its relative height com- 
pared to the bottom of the conduction band for the aver- 
age carrier density (ll|. The constant a has been chosen 
a « 0.4 in order to allow for the following approxima- 
tion: V(x,y) w for {y — yt) 2 + (x - Xi) 2 > Ar t , i.e. 
the impurity potential acts as a local perturbation. 

The effect of the scattering potential is calculated in 
the following way: first, the unperturbated wave function 
is calculated for a line x = Xi + Ari behind the impu- 
rity potential. Then the transmission probability of this 
scatterer is evaluated and the resulting wave function at 
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Ar, is calculated from: 



* T (x 4 + An, y) = $( Xi + A n , y) t{y), 



(4) 



where ^(xi + Ari,y) is the wave function, calculated 
omitting the impurity. As long as the size of the scatter- 
ing potential is much smaller than the distance between 
injector and detector QPC, the transmission amplitude 
can be calculated from: 



Kv) 
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where T is the gamma function, k 



V2m'E . 



V(xi,y) 



is the height of the impurity at the position y [Eq. (H)]; 
A = yj Arf - (yi - y) 2 (where \y - y t \ < An) is the 
half of the extension of the impurity in ^-direction for a 
given y-coordinate. Multiple impurity potentials are con- 
sidered in the same successive way. The modified wave 
function is then used to calculate the wave function at the 
detector position [x = L) using Eq. ([!]). The transmis- 
sion coefficient for the detector QPC is calculated from 
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where ^D{L,y') is the wave function in the detector 
QPC, which is taken in analogy to the wave function 
* b at the exit of the injector QPC 

From comparing the calculated transmission probabil- 
ity with the measured beam profile it is possible to de- 
termine the parameters of the scattering potential. As 
an example two experimental traces and matching theo- 
retical results are presented in Fig. ||. Good agreement 
is obtained for both cases when two scattering potentials 
are introduced into the region between injector and de- 
tector. The sizes of the impurities are comparable with 
Fermi wavelength w 0.05 /im, and their heights are 
equal or exceed Fermi energy Ep rs 8 meV. 

The sensitivity of the interference pattern to the pa- 
rameters of the scattering potential is demonstrated in 
Fig. |^ for the case of a single impurity. The initial 
parameters (location x t = 2 /mi, yi — 0.205 /im, size 
An = 0.035 /jm and height Vq ~ 16 meV) were those 
which give the best fit of the experimental data presented 



in Fig. U (curve E). In Fig. || a) the influence of the pa- 
rameter Xi on the interference pattern is demonstrated: 
if the impurity is located far away from the injector QPC 
(and detector), e.g. xt — L/2, the observed interference 
patterns are relative insensitive to small changes in Xi . A 
change of about 20% is needed to give the result shown 
in Fig. [5] a) trace 2. However, if the impurity is located 
much closer either to injector or detector [Fig. || a) trace 
3], changes in xi of 5% result in a similar strong effect 
(trace 4) as observed for the center position. The in- 
fluence of the parameter ?/, is strong but only weakly 
dependent on the x-position: changes of j/j by 5% lead 
to shifts of the interference maxima by 1 — 6% [traces 
2 and 3 in Fig. || b)]. Fig. || c) shows the sensitivity of 
the interference patterns to the variation of the param- 
eter Ari by 10%. The interference maxima for curves 2 
and 3 are shifted relative to the initial curve by 1 — 4%, 
which demonstrates also here the relatively high sensi- 
tivity of the fitting procedure. The influence of the pa- 
rameter Vo is quite different depending on the absolute 
height compared with the Fermi energy. Two examples 
are given in Fig. |3| d) where the interference patterns for 
Vq w 16 meV (trace 1) and Vo ~ 7.3 meV (trace 3) are 
shown. The lower potential height corresponds closely 
to the Fermi energy. It can be seen that a change of Vo 
of 25% (trace 2) in the first case (Vq > Ep) is need to 
achieve a similar deviation from the initial curve as it has 
been obtained for a change of only 10% (trace 4) in the 
case when Vq is comparable to Ep. 

It should be noted that the simulations do exhibit an 
ambiguity in the determination of the impurity position 
in the direction of electron beam propagation (i.e. in 
the determination of the position x{). Because of the 
symmetry of the device, a similar transmission probabil- 
ity results for an impurity located at a distance Xi or 
x^ ~ L X{ from the injector QPC boundaries, when the 
other impurity parameters are kept constant. A further 
ambiguity is in the determination of the character of the 
scattering potential, i.e., it is possible to obtain similar 
interference patterns for a barrier (Vq > 0) and for a dip 
(Vq < 0). In Fig. ^ the calculated transmission probabil- 
ity for a potential barrier (curve Bl) and a dip (curve D2) 
are shown. The results reproduce the characteristic fea- 
tures of the experimental data quite well. The scatterers 
have in both cases a radius Ar^ = 0.035 /tm but they dif- 
fer in position and height. Optimal locations found in the 
simulation were Xi — 2 /im, yi = 0.205 /im for the poten- 
tial barrier and xt = 2 /im, = 0.24 /tm for the dip, for 
potential heights are Vo ~ 16 meV and Vo w — 14 meV, 
respectively. It can be seen that a simple change of the 
sign of the potential height does not result in satisfying 
fitting of the experimental data. Fig. ^ shows the inver- 
sion from a barrier (curve Bl) into a dip (curve Dl) and 
from a dip (curve D2) into a barrier (curve B2). From 
these results it can be seen, that the absolute value of 
this fitting parameter depends mainly on the choice of 
the character of the scattering potential (dip or barrier). 
However, for most of our calculations we decided on bar- 
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FIG. 4: Calculated interference patterns for a barrier (curves 
Bl and B2) and for a dip (curves Dl and D2) in compari- 
son with experimental data (curve E). All characteristics are 
normalized with respect to the maximum transmission value. 
Results are displayed with a constant shift for clarity. Pa- 
rameters for the impurities are: x% = 2 /im, yi = 0.205 /im, 
Ari — 0.035 /mi, Vb ~ 16 meV (curve Bl); Xi = 2 /im, 
yi = 0.205 nm, An = 0.035 /im, V « -16 meV (curve Dl); 
Xi = 2 /im, J/; = 0.24 /im, Ari = 0.035 /im, Vo ~ 14 meV 
(curve B2); xi = 2 /im, j/i = 0.24 /im, Art = 0.035 /tm, 
Vo « -14 meV (curve D2). 



calculated as: 



FIG. 3: Calculated interference patterns for different values of 
the impurity potential parameters. Curve 1 in all for figures: 
Xi = 2 /im, yi — 0.205 /mi, Ar^ = 0.035 /tm, Vo » 16 meV; 
a) Xi — 2.4 /im (curve 2), Xj = 0.6 /im (curve 3) and Xi = 
0.57 /tm (curve 4); b) i/; = 0.195 /im (curve 2) and yi = 
0.215 /im (curve 3); c) Ari = 0.031 /im (curve 2) and Ar^ = 
0.039 /im (curve 3); d) Vo « 12 meV (curve 2), Vo ~ 7.3 meV 
(curve 3) and Vo ps 6.5 meV (curve 4). 



rier type potentials which correspond to 2DEG areas of 
reduced carrier density. In real GaAs 2DEG structures 
this type of scattering potentials is expected due to the 
formation of negatively charged DX centers in the remote 
donor layer 

So far we have not taken into account any temperature 
effects. Experimentally the interference patterns have 
been observed to smear out with increasing temperature, 
practically disappearing at T ~ 4 K [Fig. || a)]. We 
first consider thermal broadening as the origin of this 
behavior. The detector signal results from electrons with 
different energies Ep ± ksT propagating from injector to 
detector. In order to include this effect the calculations of 
the wave function [Eqs. (|], |[^)] were carried out for the 

j2m*(E F ±k B T) 

corresponding values ot wave vector k = — c . 

The thermally averaged transmission probability is then 



T 



df(E) 
dE 



W/2 



** D (L,y')*(L,y',E) dy' 



-W/2 



dE 



(9) 

where — is the derivative of the Fermi function jll| . 

From Fig. b) it can be seen that the reduction of the 
amplitude of the detector signal with increasing tempera- 
ture is reasonably reproduced by the thermal broadening. 
However, the signal still shows all features corresponding 
to the low temperature interference patterns [Fig. [5] b) 
solid lines], while experimentally these features are all 
smeared out. This result implies that further tempera- 
ture effects are responsible for the experimental obser- 
vations. Another phenomenon which strongly influences 
the appearance of interference effects in an electron sys- 
tem is electron-electron scattering. It leads to a dephas- 
ing of electrons but leaves their trajectories due to phase 
space restrains in a 2D system practically unaffected jl2) . 
As long as the electron-electron scattering mean free path 
l ee is larger than the distance between QPCs L we can 
use a one-collision approximation to consider this effect 
jUI, i.e. an electron that has been scattered is unlikely to 
be scattered again before it reaches the detector. It has 
been shown fL2| that the contribution of electrons that 
reach the detector ballistically can be approximated by: 
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FIG. 5: Measured electron beam profiles (a) and calculated 
interference patterns (b) at different temperatures (solid lines: 
only thermal broadening is considered, dashed lines: ther- 
mal broadening and electron-electron scattering effects are 
included). Results are displayed with a constant shift for clar- 
ity. Parameters for the scattering potentials are: Xi=l .2 ^m, 
j/i = —0.26 /im, Art — 0.08 /im, Vb ~ —8 meV for first im- 
purity (dip) and x, = 1.2 /im, yt = 0.23 fjm, Art = 0.04 ^m, 
Vo » 82 meV for second (barrier). 



$b(L,y) = #(L,y)exp ( - — 



(10) 



where l ee is determined using the usual expression for the 
energy relaxation length in a 2DEG fl6|] : 



lee = VFT e . 
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where vf — y ^§f- is the Fermi velocity, r ee is the 
electron-electron scattering time, A is the electron ex- 
cess energy relative to Ep, qxF is the 2D Thomas-Fermi 
screening wave vector, and kp = v/2 ™* Ep . Electrons 
that have been scattered lose their phase information, 
and reach the detector with an arbitrary phase. Their 
contribution to the transmitted signal can be expressed 
as: 



* s (L,y)=y°(L,y) 



1 — exp — - — 



(12) 



where ty°(L,y) is the wave transmitted from injector to 
detector QPC without impurities in the channel [Eq. (EJ)]. 
The total propagated signal is then determined as the 



sum of Eq. ( |i"o| ) and Eq. (O) . The results of this calcula- 
tion are presented in Fig.j^b) (dashed lines) for different 
temperatures from 1.7 K to 3.5 K. It can be seen that 
the effect of electron-electron scattering contributes to a 
further decrease and smearing of the interference pattern 
[Fig. H b), dashed lines]. However, even at 3.5 K the 
structure is still more pronounced than experimentally 
observed. We conclude that even more dephasing pro- 
cesses contribute to the experimental data, e.g., electron- 
phonon scattering processes, which have not been consid- 
ered in our model. 

Up to now only forward propagating electrons have 
been considered in the calculations. However, electrons 
that are scattered out of the beam can be scattered again 
by an additional impurity potential or the system bound- 
aries and then reach the detector. For a complete sim- 
ulation of the device behavior these back-scattered elec- 
trons have also to be considered in the calculations. Since 
back-scattered electrons travel much longer distances be- 
fore they reach the detector, the effect has to be taken 
into account only if the additional path length remains 
smaller than the electron-electron scattering length, i.e. 
smaller than the phase coherence length. For tempera- 
tures considered here, this holds only for back-scattering 
effects that occur at scattering centers located in close 
vicinity to the injector or detector, or if two scattering po- 
tentials are close to each other. If the impurity is located 
near the injector, the back-scattered wave must be calcu- 
lated as that part of the injected beam that is scattered 
from the impurity towards the injector QPC, reflected at 
the injector boundaries, and then retransmitted via the 
impurity towards the detector QPC. For this configura- 
tion, self-consistent calculations of the modified injected 
electron beam ^(0)' due to back-scattered electrons into 
the injector QPC should be used in the simulation. Wc 
have not attempted to incorporate such self-consistency 
into our model. However, a self-consistent method is not 
needed when the impurity is located near the detector. 
In this case, the back-scattered wave is that part of the 
beam that is reflected back from the channel boundaries 
near the detector towards the impurity and subsequently 
is scattered into the detector QPC. The calculated trans- 
mission probabilities without and with back-scattering 
effects for a temperature of T—1.8 K are presented in 
Fig. ^| a) for the case of an impurity located near detector 
QPC. It can be seen that for different values of Xi the in- 
terference patterns are very similar when back-scattering 
is omitted (lower curves in Fig. ^| a), and that the curves 
differ significantly when back-scattering is taken into ac- 
count (upper curves in Fig. ^ a). To illustrate this effect 
more clearly, a part of the upper trace in Fig. 6 a) is 
rcplottcd enlarged in Fig. 6 b) (solid line). It can be 
seen that if the initial position of the scattering potential 
is changed by an amount of the order of \f/A (Fermi 
wavelength Xp rj 50 nm for the samples under inves- 
tigation), structures in the electron beam signal which 
appeared as a maximum in the initial curve are rendered 
into a minimum and vice versa. The initial shape is re- 
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multiple back-scattering events have to be considered, 
which drastically reduce the number of electrons reach- 
ing the injector QPC. Therefore, from this observation 
no information on the propagating electron beam can be 
obtained. Furthermore, a branching of the electron beam 
can be excluded due to numerous electron beam collima- 
tion experiments usinga second QPC as detector for the 
actual electron beam |§, ffl, 17 [L8l n9| at distances of up 



to 4 /im. In none of these experiments 'dark' areas have 
been detected. The center part of the beam exhibits al- 
ways the highest signal intensity and observed structures 
are well explained by wave interference effects. It should 
be noted that the Frauenhofer diffraction picture is not 
sufficient to reproduce the observed inference patterns in 
electron beam experiments. A fully quantum mechanical 
model is needed, as we have demonstrated here. 



FIG. 6: a) Calculated interference patterns without and with 
consideration of the back-scattering effects. Parameters for 
the scattering potential are: Xi = 3.3 /im for solid lines and 
Xi = 3.4 /im for dashed lines, yi = 0.15 /im, Ari = 0.075 /im, 
Vo ~ 41 meV (the results including back-scattering are pre- 
sented with a shift), b) Enlarged section of the upper trace 
in Fig. a) for three different values of xv. 3.3 /im (solid line), 
3.313 /im (dashed line), 3.326 /jm (dot-dashed line). 

covered if the position is changed by ~ Xp/2. This effect, 
which obviously results for constructive and destructive 
interferences, directly manifests the wave nature of the 
electron beam. 

Our calculations enable us not only to understand elec- 
tron beam collimation experiments, but are applicable 
also to recent experiments where an electron beam in the 
vicinity of a QPC has been manipulated by a charged tip 
of an atomic force microscope |8j. The tip of the mi- 
croscope introduces a scattering potential which can be 
located anywhere in the 2DEG area. Electrons which are 
back-scattered into the injector are detected via changes 
in the injector current. These changes can be used to map 
out the current distribution behind the injector QPC. 
The effect of coherent constructive and destructive inter- 
ference of back-scattered electrons [Fig. 6 b)] is clearly 
demonstrated in these experiments. However, in our 
opinion the authors' conclusion of a branching of the elec- 
tron beam around the scattering potentials is a misinter- 
pretation. An experimental absence of current variation 
at the injector QPC does not imply that the current den- 
sity is zero at the point where the atomic force microscope 
tip is located. Especially when the charged tip is behind 
a potential fluctuation of the 2DEG layer no noticeable 
changes are expected in the injector current. Based on 
our simulation, we surmise that in this situation when 
the tip and the scatterer are located close to each other 



III. CONCLUSION 

In this paper we presented a comprehensive quantum 
mechanical model to calculate the transmission prob- 
ability of an electron beam traversing a 2DEG region 
confined between two opposite electrostatically defined 
QPCs. By including local potential fluctuations, tem- 
perature and geometric effects experimental results can 
be interpreted in great detail. The main features of the 
experimental data are represented in the calculation al- 
ready including only one or two locations for a circular 
shaped scattering potentials. Even though interpretation 
of the values given for the location is sometimes ambigu- 
ous due to the symmetry of the system, the extracted 
sizes and strengths exhibit only a few percent of uncer- 
tainty and give for the first time reliable values for the 
dimensional extents and strengths. It should be noted 
here that this symmetry is broken when back-scattering 
and a finite mean free path are fully taken into account. 
The extracted dimensions and amplitudes of the scatter- 
ing centers correspond well with estimations from first 
principle calculations on the correlation of remote donors 
(size: j|, strength: jui)). Also modulation of QPC cur- 
rents in local probe experiments can be completely un- 
derstood in terms of back-scattered electrons from a scat- 
tering potential induced by the charged local probe. 
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